import pandas as pd
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np



fig =  plt.figure()
ax = Axes3D(fig)
plt.rcParams['savefig.dpi'] = 600 #图片像素
plt.rcParams['figure.dpi'] = 600 #分辨率
#横坐标
xx = np.linspace(-20,20,1000,)
#纵坐标
yy = np.linspace(-20,20,1000)

X,Y = np.meshgrid(xx,yy)
Z = 1/np.sqrt(X*X+Y*Y) * np.sin(np.sqrt(X*X+Y*Y))

ax.plot_surface(X,Y,Z,rstride=1,cstride=1,cmap='hot')

plt.show()